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We investigate thermodynamic properties of a one-dimensional S = 1/2 antiferromagnetic 
Heisenberg model coupled to a lattice distortion by a quantum Monte Carlo method. In par- 
ticular we study how spin and lattice dimerize as a function of the temperature, which gives 
a fundamental process of the spin-Peierls transition in higher dimensions. The degree of free- 
dom of the lattice is taken into account adiabatically and the thermal distribution of the lattice 
distortion is obtained by the thermal bath algorithm. We find that the dimerization develops 
as the temperature decreases and it converges to the value of the dimerization of the ground 
state at r = 0. Furthermore we find that the coupling constants of spins fluctuate quite largely 
at high temperature and there thermodynamic properties deviate from those of the uniform 
chain. Doping of non-magnetic impurities causes cut of the chain into short chains with open 
boundary. We investigate thermodynamic properties of open chains taking relaxation of the 
lattice into consideration. We find that strong bonds locate at the edges and a defect of the 
bond alternation appears in the chain with an odd number of sites, which causes enhancement 
of a staggered magnetic order. We find a spreaded staggered structure which indicates that the 
defect moves diffusively in the chain even at very low temperature. 

KEYWORDS: spin-Peierls system, lattice fluctuation, quantum Monte Carlo method, bond alternation, magnetic 
susceptibility, impurity 



§1. Introduction 

The one-dimensional S = 1/2 antiferromagnetic 
Heisenberg chain is unstable to the static lattice dimer- 
ization which is well known as the spin-Peierls instability. 
The nature of the spin-Peierls transition has been inves- 
tigated extensively. Ejxperimentally apvcral compounds 
such as TTF-CuBDTE|) and CuGeOaB have been found 
out to show the spin-Peierls transition. In particular 
since the discovery of the inorganic spin-Peierls com- 
pound CuGeOs, extensive experiments have been pet- 
formed such as the neutron scattering measurementcP 
and the investigation of the effect of substitution of im- 
puritieSjO* which is difficult in the organic compounds. It 
has been found that spin-Peierls order is suppressed and 
antiferromagnetic order is favored with doping impuri- 
ties.LP Several phases appear according to the concentra- 
tion of the impurity and coexistence of spin-jfjeierls and 
antiferromagnetic states has been observed.BEl' On the 
other hand, theoretically the adiabatic lattice distortion 
has been considered and the property of the groimdialate 
has been investigated by a bosonization theory.U'H'E£F It 
is known that an arbitrarily small spin-phonon coupling 
causes spontaneous lattice dimerization in the ground 
state because the magnetic energy gain of the spin in- 
teraction in the dimerized state, which is proportional 
to u^/^ , is larger than the restoring elastic energy oc u^, 
where u is an amount of the dimerization. However, it is 
difficult to estimate thermodynamic properties exactly 
at ffirite temperature. The thermal fluctuation of the 



lattice distortion has been treated self-consistently by a 
mean field theocjijllj a Monte Carlo methocO and a 
DMRG method.la) In a mean field picture, the lattice 
has the uniform configuration above the transition tem- 
perature and the system has the same property as that 
of the uniform Heisenberg system. However, the effect 
of the thermal fluctuation of the lattice distortion is not 
clear. The quantum lattice fluctuation is considered to 
cause different effect from that of the adiabatic lattice 
distortion. It has been predicted that in case of a finite 
frequency a critical spiixL-.p|inTjnT] coupling exists for the 
spin-Peierls transition.li3'll3't£lliZP Thermodynamic prop- 
erties have been investigated by a quantum Monte Carlo 
method taking the thermaLfhictuation of the quantum 
phonon into consideration Jla^tj) It has been pointed out 
that the thermal fluctuation of the quantum phonon 
causes different thermodynamic properties in the spin- 
Peierls system from those in the static uniform chain 
even at high temperature. 

The effect of the inhomogeneity of the lattice has been 
studied extensively in quantum spin systems. It has been 
found that peculiar magnetic structures appear in the 
open chain due to the edge effect.E2lEJ For example in 
the S = 1/2 uniform Heisenberg antiferromagnet with 
open boundary, a spreaded staggered magnetic structure 
appears in the odd chaiiiyjvhile a one-node structure ap- 
pears in the even chain.l23) In the bond-alternating lat- 
tice, the defect of the alternation caaiscs a localized mag- 
netic structure around the defect £j) In the spin-Peierls 
system the change from the uniform chain to the bond- 
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alternating chain occurs spontaneously. Thus it is an 
interesting problem to study the magnetic structure in 
the open chain of the spin-Peierls system with the lattice 
relaxation. The effect of the lattice relaxation in th|jxincii 
chain has been studied in a self-consistent way.cJ'Ej'ES' 
There strong bonds locate at the edges. Furthermore, 
a localized structure of a defect of the bond alternation 
appears in the chain with an odd number of sites and a 
local staggered magnetization appears around it. Impu- 
rity problem has also been studied by the bosonization 
theoryM 

In this paper we study the effect of the thermal fluctu- 
ation of the adiabatic lattice distortion. In order to in- 
vestigate thermodynamic properties, we use a quantum 
Monte Carlo loopL-algorithm with continuous imaginary 
time (LCQMC).EIB@) This algorithm is useful to ex- 
plore the property at low temperature. Besides, there is 
no need to extrapolate with respect to the Trotter num- 
ber. We extend the standard LCQMC to the system 
with the lattice distortion. The thermal distribution of 
the lattice distortion is obtained by the thermal bath al- 
gorithm. In §2 model and method are explained. As a 
reference to the study of the thermal fluctuation of the 
lattice distortion, we review the property of the dimer- 
ized ground state investigated by an exact diagonaliza- 
tion method (ED). In §3 we investigate thermodynamic 
properties of the system with the periodic boundary. We 
show how the dimerization develops as a function of the 
temperature. In §4 we investigate properties of open 
chains, which is relevant to the effect of doping of non- 
magnetic impurities because it causes cut of the chain 
into short chains with open boundary. In §5 summary 
and discussion arc given. 

§2. Model and Method 

The one-dimensional S = 1/2 antiferromagnetic 
Heisenberg model coupled to the adiabatic lattice dis- 
tortion is given by 

N I N 

i/ = ^J.(A.)S..5.+i + -^A^ (1) 

i=l 1=1 

Ai = - u^+i, (2) 

where Ui is the lattice displacement of the i-th site from 
its equilibrium position and A^ is the lattice distortion 
owing to the lattice displacement. The distortion of the 
exchange coupling is assumed to be proportional to the 
lattice distortion 

J,(A,) = J(l + aA,). (3) 

Hereafter, for the simplicity of notation, we use the dis- 
tortion of the interaction aA^ as a variable, namely, 
Ai ~ aAi and also scaled constant k = k/a^, which 
leads 

H = ,7^(1 + A,)S, ■ + 2 E (4) 

i=l 2=1 

When we consider the pure system, the periodic bound- 
ary condition is adopted in §3. In order to investigate 
the impurity effect, we consider open chains in §4. We 



set the uniform exchange coupling J = 1 and take it as 
the unit of the energy. 

We investigate thermodynamic properties of this 
model by a Monte Carlo method. The degree of free- 
dom of the lattice is taken into account and the thermal 
distribution of the lattice distortion is obtained. The 
procedure to update the system is as follows. We start 
with an arbitrary spin and bond configuration. We up- 
date the spin configuration with the fixed bond config- 
uration by the standard LCQMC. Then we update the 
bond configuration with the fixed spin configuration by 
the thermal bath algorithm. For convenience of the sim- 
ulation, change of the bond distortion A^ is discretized 
by a small unit Aunit, which is an approximation of con- 
tinuous change. At each update of A^, three states are 
examined as the next value as depicted in Fig. 1. We 
consider a motion of the i-th site Ui. The site moves 
right or left with the unit or stays the present position, 
that is, Sui = ±Aunit or 0. According to this site move- 
ment, the bond distortions on both sides of the site are 
changed as 

J Aj_i Ai-i - 5ui 
\ A, A, + 5u^. 

In this case the total length of the lattice does not 
change: 

N 

5^A, = 0. (6) 

i=l 

In realistic materials, the total length of the lattice 
changes with the temperature and also with the mag- 
netic field. These effects are studied as the magneto- 
elastic property or the Invar problem, which are inter- 
esting topics in magnetism. However, in this paper we 
concentrate ourselves on the effect of dimering due to 
the alternate deformation of the lattice as an essential 
mechanism of the spin-Peierls phenomena. If we allow 
the total length change, the system shows the dimer- 
ization but we found that the total length decreases in 
order to gain the magnetic energy, i.e., to increase the 
amplitude of the interaction J. Thus we introduce the 
restriction eq. (||) here. The detail study of the effect 
of this restriction will be reported elsewhere. In case of 
the open chain we fix the edge sites in order to keep the 
total length of the lattice. These two updates, namely, 
LCQMC for the spin configuration and the bond-update, 
are applied alternately. 

In the simulation we set the unit of the bond distor- 
tion Aunit — 0.02 and provide the cut off A^/Aunit — 
—30 ~ 30 to restrict the exchange coupling to be anti- 
ferromagnetic {Ji = 0.4 ~ 1.6). Starting from spin all 
up configuration and uniform bond configuration, initial 
10^ Monte Carlo steps (MCS) are discarded to obtain 
thermal equilibrium state and then in most cases 10^ 
MCS are performed to sample the data of physical prop- 
erties. In case of the open chain 10^ MCS are performed 
to obtain good convergence of the data. The simulation 
through the MCS is divided into 10 bins and the errorbar 
is estimated from the standard deviation of the data of 
the bins. 
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2. 1 Dimerized ground state 

Before considering the thermal fluctuation of the lat- 
tice distortion, it would be useful to discuss ground state 
properties of the dimerized chain with 

= (7) 

in order to investigate ground state properties of the 
spin-Peierls system. Let us check the well-known relation 
So (X k~^/^. The dimerization amplitude in the ground 
state (5o is determined by ED for values of 5. In Fig. 2 
we show the 5 dependence of the change of the ground 
state energy, A£'g((5), for N = 20. There is a value of 
5 which gives a minimum value of AiJg, and the dimer- 
ized ground state is realized. In Fig. 3 we show the k 
dependence of 6o for N = 20. There we find the rela- 
tion 5q oc clearly. It is a noteworthy fact that 5q 
is afi^ectcd by the finite size effect for large k. Actually, 
(5o is found to be zero for large k in case of a small sys- 
tem size, e.g., 5q is zero for about k > 3.5 in the case 
of = 20. We have checked that the values of for 
k = 1.0 and 2.0, which will be used in the next section, 
do not change among the cases of N = 16, 20 and 24. 
Thus we confirmed the previously known properties of 
the distortion in the ground state. 

As shown in Fig. 2, it should be noted that the amount 
of the energy gain due to the dimerization AE{k) is 
rather small and also that the dependence of AEg on 
6, namely, the effective potential for the distortion, is 
very weak in the region where Aii'g < 0. Thus even at 
low temperature, e.g., at T = 0.1 for A: = 1.0, large fluc- 
tuation of the lattice is expected, which will be seen in 
the next section. 

§3. Thermodynamic Properties 

3. 1 Bond flucf,uation 

In this subsection we study the thermal fluctuation 
of the bond distortion {Aj}. First we investigate the 
distribution function of the bond distortion. In case of 
the bond alternation the distribution has two peaks at 
the corresponding position. In Figs. 4 we show the nor- 
malized distribution of the bond distortion for several 
values of k at several temperatures. As the temperature 
decreases the wide distribution becomes narrow and the 
distribution changes over from a one-peak distribution to 
a two-peak distribution. The position of the peak corre- 
sponding to the bond alternation is consistent with that 
in the ground state estimated by ED in the previous sec- 
tion. It should be noted that the distribution becomes 
asymmetrical in the intermediate stage to the two-peak 
distribution. The amplitude at the peak of positivc-A 
is larger than that at the peak of negative-A. This be- 
havior is interpreted as follows. It is easy to form the 
spin singlet state at the strong bond with positive-A so 
that the energy is lowered at the strong bond. On the 
other hand, if the spin singlet state is formed at the weak 
bond with negative-A, the energy is higher than that 
of the spin singlet pair at the strong bond. Then the 
strong bond is more stable than the weak bond. When 
the thermal fluctuation of the bond distortion is taken 
into consideration, for negative-A the effect of the ther- 



mal fluctuation is more remarkable and the distribution 

extends to wider range than positive-A. This difference 
between positive- and negative-A should disappear in 
the ground state at T = 0. In Fig. 5 we show the tem- 
perature dependence of the peak positions of positive- 
and negative-A for k = 1.0. For negative-A the peak 
appears at a lower temperature than that of positive-A 
as we mentioned above. Peak positions of positive- and 
negative-A reach to saturated values at almost the same 
temperature. 

We find that the distribution extends to wide range 
even at a modestly low temperature T = 0.1 for k = 2.0 
in Fig. 4(b). At a high temperature T = 0.5 the dis- 
tribution extends up to the cutoff (A = ±0.6). As 
shown in Fig. 4(a), the distribution extends to wider 
range for the case of fc = 1.0. The cutoff of the dis- 
tortion restricts the exchange coupling to the provided 
range (Ji — 0.4 ~ 1-6), which would represent a non- 
linear effect. If the cutoff is removed, it is expected that 
the distribution extends to ferromagnetic range and it is 
not realistic in the spin-Peierls system. Practically the 
non-linearity of the lattice distortion is expected to exist. 
Namely, the linear relation between the lattice distortion 
and the shift of the exchange coupling Jj = J(l -I- A^) 
should be changed to include higher orders for large 
A, or the non-linearity of the elastic energy potential 
U{x) = |a;^ -|- 7X^ + • • •, which implies that unreason- 
able large distortion should be prohibited. Because the 
distortion distributes almost uniformly around A = 0, it 
is expected that the value of the provided cutoff is not 
essential to represent the property of the system quali- 
tatively. Hereafter, we study thermodynamic properties 
for the case of k = 1.0, where the system shows char- 
acteristic properties of the spin-Peierls system at appro- 
priate low temperature. It is expected that thermody- 
namic properties such as the magnetic susceptibility are 
affected by the thermal fluctuation of the bond distortion 
at high temperature. 

In order to investigate the property of the bond order- 
ing, we introduce the order parameter 

which represents how the bond alternates in the whole 
chain. From the size dependence of A^g, we can estimate 
the correlation length of the alternate distortion. 

3 

- ^an (9) 

because AqAj decays exponentially and the summation 
in the above equation converges to a correlation length 
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£,{T). Thus we reach the relation 

AUN') N ■ 



(10) 



In Fig. 6 we show the temperature dependence of A^g. In 
order to investigate the size dependence, the values for 
N — 32, 64 and 128 are compared in the figure. We find 
the relation (|lO|) at high temperatures T > 0.1, which 
indicates that the bond correlation is only short range, 
^ <C A^. On the other hand, as the temperature de- 
creases, the alternate bond distortion develops near the 
temperature T ~ 0.1 where the distribution becomes to 
change into the two-peak distribution. The value of A^g 
of each size converges to the same value in the ground 
state, which indicates that the correlation extends to the 
whole chain, ^ N. In Figs. 7 we show explicitly the 
bond correlation function between the first and the j-th 
bond 



Cbond(l,i) = AiA 



J' 



(11) 



at several temperatures. We find that each bond fluc- 
tuates independently at a high temperature T — 0.5 as 
shown in Fig. 7(a). As a result, the bond distortion is 
uniform on average in the whole chain. As the temper- 
ature decreases the alternate bond correlation develops 
gradually. As shown in Fig. 7(b), at a modestly low tem- 
perature r = 0.1 short range correlation of the alternate 
distortion exists. As shown in Fig. 7(c), at a very low 
temperature T = 0.01 the alternate bond correlation is 
spread over the whole chain and the dimerized ground 
state is realized. It should be noted that the distribution 
still fluctuates from the complete dimerization. 

3.2 Magnetic susceptibility 

The magnetic susceptibility shows a characteristic fea- 
ture due to the spin gap in the spin-Peierls system. Below 
the transition temperature the energy gap opens between 
the singlet ground state and the triplet excited state as 
a result of the bond alternation. There the susceptibility 
drops exponentially with decrease of the temperature. In 
Fig. 8 we show the temperature dependence of the mag- 
netic susceptibility. In order to investigate the finite size 
effect, the values for N = 32, 64 and 128 are compared in 
the figure. The solid line denotes the susceptibility of the 
uniform chain for N = 128. We find that the suscepti- 
bility of each size drops at almost the same temperature, 
which indicates that the energy gap causing this drop 
appears on account of the bond alternation rather than 
the finite size effect. On the other hand, wc find that the 
peak position shifts to a lower temperature and the sus- 
ceptibility is enhanced at high temperatures compared to 
the case of the uniform chain. As shown in the previous 
subsection, the bond distortion fluctuates quite largely 
at high temperature. The susceptibility is also affected 
by the thermal fluctuation of the bond distortion. 

It has been pointed out that the susceptibility of the 
spin-Peierls system is different from that of the uniform 
chain at high temperature, which is due to the fact that 
the uniform phonon displacement causes the effective 
exchange coupling larger than the bare exchange cou- 



pling.t§ll3 In these studies the distortion of each bond 
behaves independently, and the total length of the chain 
is subjected to change. On the contrary, in the present 
case we fix the total length of the chain. Within this 
restriction we still found the change of the susceptibility 
at high temperature. 

In order to study the effect of the bond distortion in 
the present model on the susceptibility at high temper- 
ature, let us take into account the degree of freedom of 
the bond. Namely, we consider the following partition 
function 



N 



Z = Tr 



n 



""it) 



(12) 



where H is given by eq. (^) and the restriction on the 
fixed total length of the chain is taking into account by 
the delta-function. We approximate the density matrix 
by the Suzuki- Trotter decomposition with the Trotter 
number m = 1 



-/3J^^(l+A.)S',-S',^ 



J|g-;3J(l+A.)5..S.+ ,^ (13) 



which would be allowed at high temperature. We also 
express the delta-function by the integral formula 



2tt 



(14) 



Here we can integrate out the variables {A;} individu- 
ally 

dq 
2tt 



Tr 



n 



dA,, 



= Tr 



Xe-^'3[j(l+A,)S'.-S'.+i + |A?]+i9A. 



(15) 



where the cutoff of {Ai} is neglected in order to carry 
out the Gauss integration. Applying the Suzuki- Trotter 
formula with m = 1 again and carrying out the inte- 
gration with respect to g, we have the expression of the 
partition function 



Zoff = C'Tre"'^^"", 
/ /2\ ^ 

^ ^ i=l 

12 f N 



(16) 



2kN 



i+l 



(17) 



In this expression we find that the effective coupling 
shifts as J — > J + J^/ik, which corresponds that the 
effective temperature shifts as T T/{1 + J/Ak). This 
decrease of the temperature is inconsistent with the 
change of the peak position in Fig. 8. The second term 
Wn (Si '^i ■ '^i+i)'^ ™-^y cause shift of the peak position 
and enhancement of the susceptibility. In order to check 
these points, we investigate the susceptibility of the effec- 
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tive Hamiltonian (17) by ED. In Fig. 9 the temperature 
dependence of the susceptibiUty of this model is shown, 
which reproduces above points qualitatively. Thus we 
conclude that the model of eq. (17) represents well the 
present model at high temperature. 

As shown in Fig. 8, we find that the susceptibilities of 
the uniform chain and the present model (|^) agree with 
each other at fairly high temperature [T/ J > 0.8). How- 
ever, we find a small disagreement at high temperatures 
in Fig. 9, although the model of eq. (17) should be the 
better model in the higher temperaure. This disagree- 
ment can be understood as the effect of the restriction on 
the range of . In the simulation we restrict the range of 
the distortion A.^ , but here we allow it to change without 
restriction. Thus the susceptibility at high temperature 
depends on the details of the non-linearity of the lattice. 

§4. Inhomogeneity and Lattice Relaxation 

So far we have considered the periodic boundary con- 
dition. We found there that the lattice distortion occurs 
alternately. In the uniform Hcisenbcrg model it has been 
pointed out that the inhopLogfpeity causes various pecu- 
liar magnetic structures.E3Eil'E3) Here we consider the 
effect of the lattice distortion in open chains. It is nat- 
urally expected that the bond configuration changes in 
order to minimize the total energy of the spin system 
and the bond system. Such lattice relaxation has been 
studied by a kind of method wJifim the lattice distortion 
is determined self-consistently.oE§E3^ There a localized 
structure of a defect of the bond alternation is obtained 
and a local staggered magnetization appears around it. 
In this section we investigate how the structures of bond 
and spin develope at low temperature by means of the 
same QMC method as that in the previous section. 

In the even chain, we find that strong bonds locate 
at both ends of the lattice, and the bond alternates reg- 
ularly as shown in Fig. 10(a). Moreover, no magnetic 
structure appears. In the odd chain, strong bonds again 
locate at both ends, which causes a defect of alternation 
of the bond configuration. Namely, in lattices consist- 
ing of 4m -|- 1 sites, there appears a configuration with 
two successive weak bonds, while in lattices with 4m + 3 
sites, a configuration with two successive strong bonds 
appears, as shown in Figs. 11. In these lattices with 
an odd number of sites, it is known that a local ptnag- 
netic structure appears around the inhomogeneity.E3' In 
Fig. 10(b) we show the bond configuration for iV = 31 
at a very low temperature T = 0.01. There the profile is 
close to cosx, but not the soliton-like shape. This sinu- 
soidal shape is explaiaed by the diffusion of the localized 
soliton-like shape H' £3 Thus we conclude that the defect 
is quite diffusive and we can not determine the position 
very locally. This effect causes the average profile of the 
magnetization to be also sinusoidal as shown in Fig. 12, 
where the local susceptibility Xi is plotted. 



bility Xi°''''' defined as 
d 



h=0 



(18) 



X 



local 



dr(5f(r)5f(0)), (19) 



ss the degree of quantum fluctuation at 
^ We do not see any localized structure. 



This quantity is proportional to the local magnetization 
profile. In Fig. 13 we also show the local field suscepti- 



which indica 
each site.Ht 
either. 

Here we find that the soliton-like local defect is quite 
diflfusive and the staggered magnetic structure spreads 
in the open chain, which may cause enhancement of 
staggered magnetic order even at very low concentra- 
tion of non-magnetic doping in the spin-Peierls material, 
namely, very low concentration of the defect. 

§5. Summary 

In this paper we investigated thermodynamic prop- 
erties of a one-dimensional S — 1/2 antiferromagnetic 
Heisenberg model coupled to a lattice distortion by a 
Monte Carlo method. In particular we investigated how 
the dimerization develops as a function of the tempera- 
ture. The thermal fluctuation of the bond distortion was 
taken into consideration. 

By means of calculating the bond correlation func- 
tion, we found clear evidence that the alternate bond 
correlation develops as the temperature decreases and 
the dimerized ground state is realized at low tempera- 
ture. On the other hand, the distortion of each bond 
fluctuates independently at high temperature, where the 
alternate correlation is only short range and the distor- 
tion is uniform on average in the whole chain. Besides, 
we investigated the distribution function of the bond dis- 
tortion. As the temperature decreases the distribution 
changes over from a one-peak distribution to a two-peak 
distribution. The position of the peak corresponding to 
the bond alternation is consistent with that of the ground 
state. The distribution is asymmetric at intermediate 
temperature along the way to the two-peak distribution. 
This asymmetrical behavior is interpreted as a quantum 
effect that the spin singlet state is easy to form at the 
strong bond to lower the energy. We found that the ther- 
mal fluctuation of the exchange coupling is quite large 
at high temperature. In the present case we provided 
the cutoff of the bond distortion so that the exchange 
coupling should be antiferromagnetic. The distribution 
would extend to ferromagnetic range if the cutoff is re- 
moved. We considered that the cutoff represents a kind 
of the non-linearity, namely, the non-linearity in the re- 
lation between the bond distortion and the shift of the 
exchange coupling, or the non-linearity of the elastic en- 
ergy. Although appearance of ferromagnetic bonds is not 
realistic in the spin-Peierls system, properties of such a 
system are also of interest in a general point of view. 

We also investigated the temperature dependence of 
the magnetic susceptibility. We found that the suscep- 
tibility drops with decrease of the temperature at low 
temperature, which is a characteristic feature in the spin- 
Peierls system. We also found that the susceptibility 
is enhanced at high temperature compared to that of 
the uniform chain. In particular we found the change of 
the susceptibility even in the case where we fix the total 
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length of the lattice. 

We considered the inhomogeneity of the lattice tak- 
ing the lattice relaxation into consideration. We inves- 
tigated the structures of bond and spin in open chains. 
We found that strong bonds locate at the edges both in 
the even and odd chains, and a defect of the bond alter- 
nation appears in the odd chain. The soliton-like defect 
moves quite diffusively and the bond structure forms a 
sinusoidal configuration rather than the localized soliton- 
like one, on the contrary a soliton-like structure is ob- 
tained withiii_a_aelf-consistent treatment of the lattice 
relaxation.oSEj^ An staggered magnetic structure also 
spreads over the chain, which would give an explanation 
why very low concentration of doping of non-magnetic 
impurities causes enhancement of the staggered magnetic 
order in the spin-Peierls material.d-' 
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Fig. 4. The distribution function of the bond distortion obtained 
by QMC. (a) k = 1.0 and (b) k = 2.0 for = 128. 



Fig. 5. The temperature dependence of the peak position of the 
distribution of the bond distortion obtained by QMC. fc = 1.0 
and TV = 128. The symbols open circle and solid circle denote the 
peak position in positive and negative range of A, respectively. 
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Fig. 7. The bond correlation function C^,oiid(^T j) = ^l^j ob- 
tained by QMC. (a) T = 0.5, (b) T = 0.1 and (c) T = 0.01 for 
k = 1.0 and N = 128. 



Fig. 8. Tlie temperature dependence of tlie magnetic susceptibil- 
ity obtained by QMC. k = 1.0. Tlie solid line denotes the values 
of the uniform Heisenberg chain for N = 128. 



Fig. 9. The temperature dependence of the magnetic suscepti- 
bility of the effective Hamiltonian (17) obtained by ED. k = 1.0 
and N = 12. The solid line denotes the values of the uniform 
chain. 



Fig. 10. The bond configuration in open chains obtained by 
QMC. (a) Af = 32 and (b) Af = 31 for fc = 1.0 at T = 0.01. 



Fig. 11. The conceptual configuration of the bond alternation 
with a defect in odd chains with (a) 4m -|- 1 sites and (b) 4m -|- 3 
sites. The thick line denotes the strong bond. 



Fig. 12. The local susceptibility Xi obtained by QMC. For the 
open chain with iV = 31 for fe = 1.0 at T = 0.01. 



Fig. 13. The local field susceptibility x'"'"'"' obtained by QMC. 
For the open chain with JV = 31 for fe = 1.0 at T = 0.01. 
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